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^ . 1 Introduction 

Q^ ', The experience with data analysis in inclusive measurements of deep inelastic scattering (DIS) 

^ I shows that radiative effects are very important and should be taken into account with maximally 

^ I possible accuracy. Numerical results obtained for the case of the HERA collider ||l|, |, ||, ||, ||, ||, 

^ § have shown that the ratio of the observed and the Born cross sections varies from 0.5 to 2.0, 

^' and even larger values are reached at the borders of the kinematical region. The calculation of 

r-| ! radiative corrections (RC) for semi-inclusive and exclusive reactions can not be easily reduced 

to the case of inclusive scattering, but rather requires separate considerations. There are two 

k>( ; important approaches which can be applied for this task. The first one is the so-called semi- 



a 



^ ■ analytical approach, where covariant formulae for RC are obtained after applying the procedure 
of covariant cancellation of infrared divergences and integration over the photon phase space. 
No other approximation than the ultrarelativistic one is used within this approach. Explicit 
analytical formulae for RC in diffractive vector meson electroproduction were obtained in Ref. 
]T0|. The radiative correction factor has the following form 



V = -^ = exp((5i„/)(l + S„. + 5yac) H ^^- (1) 

Here S^r and S^ac come from so-called virtual corrections: the vertex function and vacuum 
polarization by lepton pairs and hadrons, resp. 6inf appears after cancellation of infrared 
divergences. The exponential results from the summation of soft photon contributions over 
all orders of perturbation theory, ahard is the contribution of hard photon radiation. In the 
general case it has the form of a three-dimensional integral over the photon phase space. In 
Ref. [0 the numerical analysis was done for fixed target experiments. In this report we provide 



numerical results also for deep inelastic scattering at HERA. Moreover, we discuss the case of 
photoproduction within the semi-analytical approach. 

Another possible approach for the calculation of RC in semi-inclusive and exclusive processes 
is based on using a Monte-Carlo generator for inclusive DIS. In this report we discuss how the 



RC can be calculated using information obtained during event generation. Also we show that, 
under certain assumptions, both methods give the same results. This fact is illustrated for the 
case of measurements at HERMES. 

It should be noted that all numerical results are shown for the electron method of the 
reconstruction, when the kinematical variables are calculated using the measurement of the 
final lepton. There are other reconstruction methods (for example, the double angle method 
using for the case of diffractive vector meson production), which are, in the most of the phase 
space, less sensitive to these radiative corrections. 



2 The code DIFFRAD 

The FORTRAN code DIFFRAD was developed to calculate radiative corrections for diffractive 
vector meson electroproduction and is based on the semi-analytical approach. This code has two 
versions which are referred to as IDIFFRAD.F and MDIFFRAD.F. The first version allows to 
calculate the RC to the cross section [d^a / dxdQ'^dtd(f)h and d^a/dxdQ'^dt) of diffractive vector 
meson electroproduction using usual (non Monte-Carlo) methods of numerical integration. The 
second one exploits the same formulae but using Monte Carlo methods for integration. It 
allows to integrate not only over the photon kinematic variables, but also over those of the final 
lepton. Integration over Q^, W^ and t in arbitrary bins can be performed as well as over the 
full kinematically allowed phase space of the final lepton. Therefore MDIFFRAD.F has more 
possibilities in comparison to IDIFFRAD.F, but it requires more CPU time. 

The code requires input for the model of the virtual photoproduction cross sections a^^x and 
allows to apply a cut on inelasticity. As input for ai^T we use the model presented originally in 
Ref. [Til and refined in Ref. ||12[ (DIPSI). The implementation of other models is possible and 
straightforward. A cut on the inelasticity v {v = p\ — M^, where ph is the total momentum 
of unobserved particles and M is the proton mass) is often suitable to reduce the RC. The 
value of the cut is related to the experimental resolution on v and depends on details of the 
experimental procedure of non-exclusive background subtraction. A detailed discussion of this 



point can be found in Ref. ||10[- Both versions of DIFFRAD allow one to calculate RC with 
and without this cut. Note that two additional options also exist for the calculation of the 
inelasticity distribution with and without RC and for the fast evaluation of an approximate 
calculation of ahard- 

The FORTRAN code DIFFRAD can be received upon request via email from aku@hep.by 



3 Electroproduction 

Results obtained with the help of DIFFRAD for the three-fold differential cross section are 
presented in Fig. |l] (-y/i=300 GeV). Different plots in this figure show that there is no strong 
VT-dependence of the RC factor tj. This is due to the fact that ai and aT are almost fiat in 
the considered kinematical region. Moreover there is no large explicit dependence on W in the 
RC factor itself. The quantities 5^r and 5^ac in Eq- (0) are practically independent of W and 
only a small dependence comes from 5inf and ahard- The Q^- dependence shown in this figure 
is typical for the case when the inelasticity cut is not applied. If this cut is used then the 
rise of t] when Q^ goes down would be suppressed. From this figure one can also see that the 
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Figure 1: RC factor for p(770) electroproduction at HERA; y^=300 GeV. Symbols from top 
to bottom correspond to t = —0.7, —0.5, —0.3, —0.1 GeV^ 



t-dependence is rather important. Figures and ^ illustrate this last property, rj crosses unity 
for —t ~ 0.25-0.3 GeV^ and rises with increasing |t|. The large positive correction in this case 
is a result of the large phase space for photon radiation. The cut applied to the inelasticity 
(or E — pz) can again reduce the value of the RC factor for large \t\. As a consequence of the 
t-dependence of t], the observed slope parameter also receives large RC. This is illustrated in 
Fig. ^ The Born cross section has a steeper slope with respect to |t| than the observed cross 
section. The radiative correction to the slope parameter is negative and ~ 10%. 

The model used for exclusive virtual photoproduction of vector mesons is based on the 
hypothesis of s-channel helicity conservation (SCHC)[]. In the case of SCHC the Born cross 
section is independent of the azimuthal angle (ph between the scattering and the production 
plane. However, the observed cross section does depend on (ph, because the formulae for RC 
includes this angle. Figure || shows an example. Since only atard has an essential dependence 
on the azimuth, r] has a visible (ph dependence at large |t|, where atard gives a relatively large 

^We should note that now there are experimental indications for the violation of this hypothesis |l^, Q as 
well as theoretical models beyond SCHC |15, |^. Unfortunately the development of DIFFRAD for using such 
models is not straightforward, because the calculation of additional terms is required. The calculation for the 
general case of scalar meson production is considered in p7| . 
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Figure 2: Observable and Born cross 
sections for exclusive p(770) electro- 
production as a function of t; y^=300 
GeV; 1^=70 GeV; Q^=A GeV^. 



6,GeV-^ 








• 
• 


• 


' • • . 


• 


■ m 


■ 


■ - . 

■ 


■ 


: 




® - Born slope 








■ - observed slope 








< < 1 < < < < 1 < < < < 1 < < < < 1 < < 





0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

-t,GeV2 

Figure 3: Slope parameter for 

p(770) electroproduction; ^/s=300 
GeV; W=70 GeV; Q^=A GeV\ 



contribution to the total RC. 

We note finally that the RC to the cross section of diffractive vector meson electroproduction 
is very sensitive to the inelasticity cut. This cut suppresses the contribution of hard photon 
radiation, which is always positive. Thus the use of harder cuts leads to smaller values of the 
radiative correction factor. However, in the region of small |t| where radiative corrections are 
negative (r/ < 1) the influence of the cut is not essential or leads to larger absolute values of 
RC (^-1). We do not apply the cut for this numerical analysis. 

All numerical results shown here were obtained for the case of exclusive p(770) production. 
However there is no essential dependence on the type of the observed vector meson, and all 
discussed features are very similar for the production of heavier vector mesons. 



4 Photoproduction 

In photoproduction interactions, i.e. very small Q^ (quasi-real photon exchanged), the final 
positron is scattered at very small angle and escapes detection in the main detector through 
the beam pipe. In this case QED radiative corrections can be calculated by integration of the 
analytical formulae over the phase space of the final positron. A numerical analysis shows that 
radiative corrections in this case are smaller than in the case of electroproduction. The results 
are shown in Fig. ^. Note that small values of RC in this case reflect the statement of the 
Kinoshita-Lee-Nauenberg theorem [^, |19| . For our case it says that if we integrate over the 
ice of a particl 
20[] , for example). 



phase space of a particle which radiates, then all leading logarithms are mutually cancelled (see 
Ref. ^ 
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Figure 4: RC factor for exclusive 
p(770) electroproduction as a function 
of 0^; yi=300 GeV; W^=70 GeV. 



Figure 5: RC factor for exclusive 
p(770) photoproduction as a function 
of t; v^=300 GeV. 
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Figure 6: Possible scheme of Monte 
Carlo calculation of the RC factor. 
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Figure 7: RC factor for HERMES 
from DIFFRAD (semi-analytic ap- 
proach, full line) and from RADGEN 
(Monte Carlo approach, points with 
error bars); y/s=7.18 GeV; Q^=1.5 
GeV^; W^=20 GeV\ 



5 Monte Carlo approach 

The Monte Carlo approach considered here is based on a generator for inclusive processes with 



as 



radiative effects included. Here we focus on the case of HERMES and use RADGEN [^ 
such a generator. This generator was developed for the fixed target experiments while for 
HERA experiments HERACLES |2^ which includes electroweak effects should be used. 



Two Monte Carlo samples with and without radiative effects have to be prepared. The 
procedure is sketched in Fig. |^. Firstly, the kinematical variables v and Q^ defining the mo- 
mentum of the final lepton are generated in accordance with the Born cross section. In the 
case of generation of the sample without RC these variables are used as an input to generate 
the momentum of the vector meson. For the generation of the sample with RC, a Monte Carlo 
generator including real photon emission has to be used. Three main tasks are performed by 
this generator: 

1. simulation of the appropriate scattering channel (non- radiative; elastic, quasielastic or 
inelastic radiative tails) in accordance with their contribution to the total observed cross 
section; 

2. calculation of the inclusive RC factor r^*"'^ for the given values of v and Q^ . Since for both 
samples we generate the same number of events this factor has to be used as a weight to 
obtain the observed cross section as a weighted sum of events^; 

3. generation of the kinematical variables Vtme and Q^^^e which are then used as input for 
generating the momentum of the vector meson. Note that Vt^^e and Q%ue differ from z/, 
Q^ for events which contain a radiated photon. 

The events are collected in predefined bins. If the numbers of events in a given bin for the 
samples with and without radiative corrections are n\ and n2, resp., the RC factor for vector 
meson electroproduction r]p can be estimated as 

,. = ■'-<'^> (2) 

where t]*"*^ is the RC factor for the inclusive process. We note that due to photon radiation 
^tTue < ^- Thus the exclusive production of a vector meson with a given value for its energy 
(which is related to |t|) is less probable for a radiative event than for a non-radiative event with 
the same meson energy. As a result n\ < n2 and the typical value of the RC factor for vector 
meson production (for semi- inclusive processes in general) is smaller than for the inclusive case. 

Figure ^ presents results for the HERMES experiment obtained within the two approaches. 
From this picture we can conclude that the Monte Carlo approach correctly reproduces the 
RC factor for vector meson electroproduction. The inelasticity cut was not applied for these 
calculations. We note that the parameter A which separates soft from hard photon radiation 
and which is always necessary in the Monte Carlo approach (see discussion in Ref. |^1]]) should 



be chosen smaller than its default value A = lOOMeV. This value is close to the threshold energy 



^Another algorithm without calculation of weights is possible. In this case both samples with and without 
RC are generated in accordance with the observed and the Born cross section, resp. Both approaches are 
equivalent, and everything discussed below is valid for this case as well. 



for observing photons in the HERMES detector. In the case of vector meson electroproduction 
the maximal photon energy is smaller (especially for small |t|) in comparison with the inclusive 
case. So the parameter A has to be reduced also. For the present calculation A = 5MeV was 
chosen. For smaller values of A we do not see an essential dependence on A. 
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